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One of the most primitive but elusive current-voltage (I-V) responses of a superconductor is when 
its supercurrent grows steadily after a voltage is first applied. The present work employed a mea- 
surement system that could simultaneously track and correlate I{t) and V{t) with sub-nanosecond 
timing accuracy, resulting in the first clear time-domain measurement of this transient phase where 
the quantum system displays a Newtonian like response. The technique opens doors for the con- 
trolled investigation of other time dependent transport phenomena in condensed-matter systems. 
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A particle under the action of a single applied force 
accelerates ballistically in accordance with Newton's sec- 
ond law. In the presence of a frictional force, an applied 
force will ultimately maintain a constant velocity rather 
than produce acceleration. Analogously, an externally 
applied voltage V maintains a constant current / in the 
case of a resistive conductor, whereas it can ballistically 
accelerate the superfluid in a superconductor, leading to 
a supercurrent that grows with time. This acceleration 
phase of the supercurrent lasts for a very brief period — 
until flux motion, phase slip centers, or other dissipative 
phenomena set in — making it extremely difficult to ob- 
serve in the time domain in a correlated current-voltage 
measurement. In the present work, an electrical measure- 
ment system was developed that could resolve and corre- 
late the time evolutions of /(i) and V{t) on a subnanosec- 
ond time scale. For the sample pattern, extremely long 
meander geometries (with length-to-width aspect ratios 
in the thousands) were employed to prolong the acceler- 
ation time while maintaining y at a manageable level. 
The combination of these two measures facilitated the 
successful observation of the acceleration phase. 

The acceleration of the supercurrent density js is given 
by (from the London equations QjQ) 
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where E is the local internal electric field, e is the elec- 
tronic charge, m* is the effective mass and is the super- 
fluid density (related to the number of electrons per vol- 
ume participating in the condensate); the far right hand 
side of the equation relates Us to the London magnetic- 
field penetration depth A l we can take j = js + jn — js 
because the normal current density jn is a negligible com- 
ponent of the total current density j. This supercurrent 
acceleration phase lasts for the duration At « jc^x^X^/E, 
where jc is the critical current density that marks the 
onset of resistance. The inductance-like proportionality 
between dj/dt and E in Eq. [1] arising from the inertia of 
the superfluid, is referred to as the kinetic inductance Lk- 
In terms of the geometrical length I and cross sectional 



Lk = 
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where A is a more general penetration depth, which in- 
cludes effects such as impurity scattering (A > A^). 
Kinetic inductive effects are small except close to the 
transition temperature Tc, where their signatures have 
been seen in the high-frequency ac response or as non- 
equilibrium inductive voltage spikes during abrupt cur- 
rent steps [3-ll|. In the present work, timescales were 
chosen to be short enough to have a sufficient magni- 
tude of V while long enough (compared to characteristic 
timescales such the gap-relaxation and electron-phonon 
scattering times) to avoid non-equilibrium effects. Vari- 
ations in fields occured at length scales that were long 
compared with both A and the coherence length ^, so 
as to avoid non-local effects. Thus the conditions were 
optimum for observing the simplest limiting behavior of 
an accelerating condensate as predicted by the London 
equations, i.e., Eq. [TJ 

The samples used in this work were niobium films de- 
posited on silicon substrates with DC magnetron sputter- 
ing. The films were patterned into long narrow meanders 
by electron-beam lithography using the lift-off technique. 
Sample A had a thickness of i = 70 ± 8 nm, a width of 
w = 12.1±0.6 /im, and a length (between voltage probes) 
of / = 4.80 ± 0.01 cm. Sample B had the dimensions 
< = 85 ± 8 nm, u; = 8.9 ± 0.6 ^i^Jl, and I = 4.53 ± 0.01 
cm. Their respective superconducting transition temper- 
atures were Tc=6.74 K and Tc=7.23 K. 

The measurements were carried out in a pulsed-tube 
closed-cycle refrigerator in zero applied magnetic field. 
The electrical measurements were conducted with a 
pulsed signal source and detection electronics, in com- 
bination with a digital storage oscilloscope. Parts of 
the signal-source and preamplifier circuitry in this setup 
were developed and built in-house. An active (buffered) 
ground arrangement was developed for improving the 
shielding between the fast changing high-voltage signal 
in the current leads from the low-voltage sample-signal 
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sensing leads. The entire signal chain up until the digi- 
tal oscilloscope is analog. Using pulsed signals instead of 
continuous ac or dc excitations permits a wider range of 
currents without Joule heating of the sample and a flex- 
ible control over the waveform shape. The system per- 
forms simultaneous independent differential four-probe 
measurements of I{t) and V{t) with a relative timing ac- 
curacy of ^100 ps. The stray mutual inductive coupling 
between current and voltage leads has a (temperature 
independent) value of < 1 nH (the self inductances of 
the leads themselves are not sensed because of the four- 
probe configuration). The absolute accuracy of the in- 
ductance values measured in this system is about ±5%. 
The voltage-measurement sensitivity is about 1 iiY. The 
time interval between digitized samples is 10~^° second 
(the single-shot digitizing sampling rate is 10 GS/s). The 
speed and accuracy with which both I{t) and V{t) were 
tracked and correlated in a superconductor in the present 
experiment are, to our knowledge, unprecedented. 

Some tests and verifications of the measurement sys- 
tem are shown in Fig. [TJ Panel (a) shows the voltage and 
current (scaled by a constant) for a purely resistive test 
sample and panel (b) shows the current derivative and 
voltage (scaled by a constant) for a purely inductive test 
sample. The time scales used in the actual experiment 
were longer than these test conditions of Fig. [T] so that 
the temporal tracking between the current and voltage 
sensing circuits was essentially perfect. Some additional 
information on the apparatus can be found in our previ- 
ous review articles [H, • 
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FIG. 1: Measurement- apparatus temporal accuracy checks, 
(a) V{t) and I{t) (multiplied by a constant R=62 ^) for a 
test resistor in place of a superconducting sample, (h) dl{t)/dt 
and V{i) (divided by a constant L=15 nH) for a test induc- 
tor in place of a superconducting sample. The voltages across 
the purely resistive and inductive loads are seen to track their 
respective current and current-derivative functions with sub- 
nanosecond accuracy. 



Figs. [2ja) and (b) show V{t) (solid lines) across two 
niobium-meander samples in the superconducting state 
at one temperature. Panels (c) and (d) show the cor- 
responding I{t) functions, which are seen to accelerate 
steadily during the plateaus in V{t). The dashed lines 
in panels (a) and (b) show dl/dt scaled by a constant 
(i=16.7 and 16.9 nH for samples A and B respectively) 
and are seen to track V{t) in instantaneous detail. Thus 
the response is purely inductive, with an inductance that 
is independent of / and dl/dt. 
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FIG. 2: Time dependencies of voltage and current in niobium 
meanders at T=3.84 K. Panels (a) and (h) represent V{t) 
and dl/dt (scaled by a constant L) for samples A and B re- 
spectively. Panels (c) and (d) show the corresponding I(t) 
functions, which rise steadily during the voltage plateaus of 
the panels above; the time axes for panels (a) and (c) and for 
panels (b) and (d) are the same. 

The ratio between the dl{t)/dt and V{t) curves in the 
top panels of Fig. [2] gives the time and current depen- 
dent inductance: L{t) = V{t)/[dl{t)/dt]. Figs. El; a) and 
(b) plot this L[t) versus time for various values of T, for 
each of the samples. The plateau value of L is seen to in- 
crease steadily with temperature as is expected because 
of the declining superfluid density and consequent rising 
A L ■ Another interesting trend is that the curves at high- 
est temperatures show L{t) functions that rise with time 
(i.e., current). This happens because the current sup- 
presses the superfluid density through its pair-breaking 
action, a regime not seen before in any other kind of 
measurement. Note that continuous-ac probes of rig can- 
not endure high enough excitation levels to explore this 
regime because of Joule heating; and tunneling measure- 
ments reveal the spectral gap Q,g rather than n^. A sys- 
tematic study of the suppression of by j will be the 
subject of a future investigation, since the optimum sam- 
ple geometry for studying this effect is almost opposite to 
the sample geometry required for the present experiment. 

Figs, ^c) and (d) plot L (as measured above) versus 
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FIG. 3: (a) Measured total inductance, L{t) — V (t) .dt / dl (t) , 
versus ttme for sample A. The curves correspond to the tem- 
peratures (from top to bottom): T — 6.10, 5.91, 5.54, 5.17, 
4.73, 4-10, and 3.78 K. (Each plot symbol on these curves cor- 
responds to separately measured digital voltage sample. The 
period between samples is 100 ps.) (b) A similar L{t) plot 
for sample B. For this panel, the temperatures (from top to 
bottom) are: T = 6.51, 6.4, 6.22, 6.02, 5.84, 5.24, 4-64, and 
3.84 K. Panels (c) and (d) show the temperature dependen- 
cies of the above measured L values (taken on the plateaus 
around t ~ 11 ns) for samples A and B respectively. The 
symbols show the experimental data and the solid line rep- 
resents the least-squares fit to the two-parameter function 
L{T) = Lg + Lfc(0)/[1 - [T/Tcf] as discussed in the text. 



T for each sample. This total inductance L =Lk +Lg 
has components corresponding to the kinetic inductance 
Lk as well as a geometrical inductance Lg (from mag- 
netic flux change). The temperature dependence arises 
chiefly from Lk', the changes in Lg with temperature — 
arising from changes in the current-density profile across 
the cross section — are relatively small (~4%). From Eq.[5] 
and the empirical temperature dependence of the pene- 
tration depth we have 
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where ^^(0) = fioX'^{0)l/A. The solid line curves in 
Fig. ^c) and (d) correspond to a least-squares fit to the 
function L(T) = Lg + L/c(0)/[l - (T/tJ)^]. The values 
of Lk{0) and Lg obtained from this fitting {Tc is not a 
fitting parameter) are listed in columns 2 and 4 of Ta- 
ble m The coefficients of determination of the fits are 
i?2=0.9989 and i?2^o.9994 for samples A and B respec- 
tively. The standard errors of the fit combined with the 
error in the inductance measurement gives the error bars 
for Lk{0) that are indicated in the table. 

The third and fifth columns of Table H] show, for com- 
parison, theoretical estimates of Lfc(O) and Lg. For find- 
ing Lk (0) we note that the effective penetration depth (A) 
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2.8±0.2 
2.8±0.2 


3.5±0.7 
3.7±0.7 


12.2±1 
12.9±1 


15.4±3 
12.0±3 



TABLE I: Experimentally observed values and theoretically 
estimated values of the kinetic and geometrical inductances. 



becomes lengthened with respect to its clean-limit value 
(Al) in the presence of scattering by static disorder. This 
effect of impurity scattering can be incorporated through 
the residual resistivity Po and expressed in terms of the 



order parameter A as [14 1 
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Taking the measured values of /9o=0.347 /iJ^.m and 0.369 
/iil.m, and obtaining A(0) w iksTc from our measured 
values of Tc, we get A=:223 nm and 222 nm, and Lk—i.b 
nH and 3.7 nH for samples A and B respectively. There 
is an uncertainity in the values po because of the un- 
certainity in sample dimensions and an uncertainity in 
A(0) because of an uncertainity in the absolute value of 
Tc (this is roughly estimated to be around 200 mK). This 
gives rise to the error bars in the theoretical Lk values 
that are stated in the table. The calculated values of Lk 
are somewhat larger than the measured ones of but of 
comparable magnitude. 

The theoretical geometrical inductances for the mean- 
ders, tabulated in the last column of Table HI were com- 
puted numerically by integrating the magnetic flux that 
links to the path between the voltage probes. The error 
bars in the theoretical Lg arise from the uncertainities in 
the sample dimensions and the approximations inherent 
in the calculation. It can be seen that the theoretically 
estimated Lg values are also in agreement with their mea- 
sured counterparts. 

In summary, this work has explored the initial accel- 
eration phase during which a supercurrent builds up in 
response to an applied voltage. The voltage and current 
curves of Fig. [2] represent the first clear and direct time- 
domain demonstration of this primitive regime, where the 
quantum system shows a Newtonian like response. It is 
also the first time to observe the non-linear regime where 
the current suppresses the superfluid density, thereby in- 
creasing the kinetic inductance. The instrumentation 
developed for this experiment is unique and represents 
the first measurement of its kind where both V{t) and 
I{t) are tracked in a superconductor with sub- nanosecond 
timing accuracy. This technique can reveal more detailed 
information than just an impulse-response measurement, 
and it can be used to explore time-dependent and non- 
equilibrium phenomena in condensed-matter systems in 
a controlled way (some examples of such regimes in su- 
perconductors would be those related to phase slippage, 
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glassy dynamics, and the nascent stage of a vortex right 
after its nucleation). The present work and its method 
should be distinguished from past experiments in which 
an abrupt supercritical current pulse was applied [llj 
and only the subsequent V{t) response was measured 
without monitoring I{t). In those experiments the su- 
perconductor recoils in a highly non-equilibrium man- 
ner to the supercritical stimulus. In the present study, 
the superconducting system is always maintained close to 
equilibrium by keeping the experimental timescales well 
in excess of the gap-relaxation and electron-phonon scat- 
tering times, while keeping the timescales short enough 
to observe the inertia of the superfluid. 

The authors acknowledge useful discussions with J. M. 
Knight, B. 1. Ivlev, R. A. Webb, R. J. Creswick, T. R. 
Lemberger, G. Simin, T. M. Crawford, and F. T. Avi- 
gnone III. This research was supported by the U. S. 
Department of Energy through grant number DE-FG02- 
99ER45763. 



[lo; 

[11 



[12: 

[is: 



* Corresponding author email: kunc hur@sc.edu1 [14] 
URL: http : / /www . physics .sc. edu/kunchurl 
[1] F. London and H. London, Proc. Roy. Soc. (London) 



A149, 71 (1935). 

"Introduction to Superconductivity" by M. Tinkham, 
2nd Edition. McGraw Hill (New York, 1996). 
I. F. Oppenheim, S. Frota-Pessona, and M. Octavio, 
Phys. Rev. B. 25, 4495 (1982). 

A. Geier and G. Schon, J. Low Temp. Phys. 46, 151 
(1982). 

D. J. Frank, M. Tinkham, A. Davidson, and S. M. Faris, 
Phys. Rev. Lett. 50, 1611 (1983). 

S. M. Anlage, H. S. Howard, J. Snortland, S. Tahara, 

B. Langley, C. B. Eom, and M. R. Beasley, Appl. Phys. 
Lett. 54, 2710 (1989). 

S. M. Anlage, B. W. Langley, G. Deutscher, J. Halbritter, 
and M.R. Beasley, Phys. Rev. B 44, 9764 (1991). 
J. Y. Lee, and T. R. Lemberger, Appl. Phys. Lett. 62, 
2419 (1993). 

S. D. Brorson, R. Buhleier, J. O. White, I. E. Trofimov, 
H. -U. Habermeier, and J. Kuhl, Phys. Rev. B 49, 6185 
(1994). 

S. Cho, J. of the Korean Phys. Soc. 31, 337 (1997). 

F. S. Jehia, J-P. Maneval, F-R. Ladan, F. Chibane, 

A. Marie-de-Ficquelmont, L. Mechin, J-C. Villegier, M. 

Aprih, and J. Lesueur, Phys. Rev. Lett. 81, 1933 (1998). 

M. N. Kunchur, Mod. Phys. Lett. B. 9, 399 (1995). 

M. N. Kunchur, J. Phys.: Condens. Matter 16, R1183 

(2004). 

T. R. Lemberger, I. Hetel, J. W. Knepper, and F. Y. 
Yang, Phys. Rev. B 76, 094515 (2007). 



